Associations of IFT20 and GM130 protein expressions with clinicopathological features and survival of patients with lung adenocarcinoma

Background Lung cancer is the leading cause of malignancy-related mortality and lung adenocarcinoma accounts for about 40% of lung malignancies. The aim of this study was to investigate the associations of intraflagellar transport protein 20 (IFT20) and Golgi matrix protein 130 (GM130) expression with clinicopathological features and survival in patients with lung adenocarcinoma. Methods The expressions of IFT20 and GM130 protein in cancerous and matched adjacent lung tissues of 235 patients with lung adenocarcinoma were assessed by tissue microarray and immunohistochemistry, which were indicated by the mean optical density (IOD/area), the rate of positive staining cells and staining intensity score. The correlation between IFT20 and GM130 protein was assessed by Spearman’s rank correlation. Associations of IFT20 and GM130 protein expression with clinicopathological features of patients were analyzed by multivariate logistic regression models. The survival analysis of patients was performed by Cox proportional hazard regression models. Results With adjustment for multiple potential confounders, each one-point increase in IFT20 protein staining intensity score was significantly associated with 32% and 29% reduced risk for TNM stage in II ~ IV and lymphatic metastasis of patients, respectively (P < 0.05). And each one-point increase in GM130 protein staining intensity score was associated with a significant reduction in the risk of poor differentiation and tumors size > 7 cm by 29% and 38% for lung adenocarcinoma patients, respectively (P < 0.05). In stratified Cox model analysis, enhanced IFT20 staining intensity score was significantly decreased the risk of death by 16% for patients without distant metastasis. And elevated the IOD/area of GM130 expression significantly decreased the death risk of lung adenocarcinoma patients with tumor size > 7 cm or distant metastasis by 54% and 65%, respectively (P < 0.05). Conclusion IFT20 and GM130 protein expressions were negatively associated with tumor differentiated types, size, TNM stage and lymphatic metastasis of lung adenocarcinoma. Both IFT20 and GM130 proteins have some protective effects on the survival of lung adenocarcinoma patients with specific clinicopathological features. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-022-09905-6.


Introduction
Lung cancer is the leading cause of malignancy-related mortality, resulting in more deaths than breast, prostate, colorectal and brain cancers [1]. Lung adenocarcinoma is the most common histological subtype in non-small cell lung cancer (NSCLC), which accounts for about 40% of lung malignancies [2]. The average 5-year survival rate of lung adenocarcinoma is only 18%, mainly due to invasion and metastasis of tumor cells, diagnosed at advanced stages or lack of effective treatment strategies [3]. The tumor invasion and metastasis have a significant impact on the quality of life and overall survival (OS) of patients with lung adenocarcinoma [4]. Thus, early detection of tumor invasion and metastasis is crucial to improve the survival rate of patients with lung adenocarcinoma.
Primary cilia are antenna-like organelles protruding from the surface of most eukaryotic cells that participate in cell proliferation, differentiation and migration [5]. Some studies have reported that the number of primary cilia markedly decreased in the transformation of many tumor cells, including breast, pancreas, ovarian, and kidney cancer cells, whereas it was significantly higher in cancerous tissues of lung adenocarcinoma, colon adenocarcinoma and follicular lymphoma than that in normal control tissues [6]. Primary cilia require a series of intraflagellar transport (IFT) protein complexes that participate in the transportation of ciliary signaling proteins as well as the process of cellular protrusion and resorption to assemble and maintain their structures [7]. IFT20, the smallest IFT protein, localizes at the cis-Golgi apparatus that regulates the assembly and maintenance of primary cilia by sorting of ciliary cargo from the cis-Golgi to the base of the cilium [8,9]. Down-regulated expression of IFT20 mRNA or protein could not only inhibit ciliary assembly and decrease the quantity of primary cilia in mammalian cells, but also reduce invasion and metastasis of colorectal cancer and osteosarcoma cells [8][9][10]. A new study found that IFT20 promotes collective invasion of colorectal cancer by regulating organization of Golgiassociated, stabilized microtubules and Golgi polarity in colorectal cancer cells at the invasive front [9].
Recent studies indicated that tumor's invasion and metastasis often require cell polarization and migration [10], which need the Golgi to adopt a ribbon-like morphology in promoting polarized secretion to the cell surface [11,12]. As a cis-Golgi matrix protein, GM130 is essential for the maintenance of Golgi structure and protein transportation that involves in cell polarization and migration [13,14]. Expression of GM130 mRNA or protein was significantly lower in tissues of colorectal adenocarcinoma and breast cancer than that in matched normal tissues, and the depletion of GM130 mRNA enhanced the invasiveness of breast cancer cells [15,16]. Remarkably, IFT20 protein localized at Golgi could combine with GM130 and A-kinase anchor protein 450 (GM130-AKAP450) complex to regulate cellular microtubules nucleation, which contributes to Golgi ribbon formation in achieving polarized secretion for cell migration and invasion [10,13]. Thus, both IFT20 and GM130 proteins might play important roles in the invasion and metastasis of various malignant tumors, but their roles in the tumorigenesis and development of lung adenocarcinoma remains unclear so far.
Hence, we determined the expressions of IFT20 and GM130 protein in cancerous and matched adjacent lung tissues of patients with lung adenocarcinoma by tissue microarray and immunohistochemistry to investigate their potential roles in the development of lung adenocarcinoma and the relationships between their expressions and survival of patients with lung adenocarcinoma.

Study population
A total of 235 patients with lung adenocarcinoma were selected from one hospital in Wuhan city in China from March 2010 to September 2015 by simple random sampling, and followed up to June 2018. All the patients underwent lung cancer surgery and were finally diagnosed with lung adenocarcinoma by pathological biopsy. This study was approved by the Ethnics and Human Subject Committees of Tongji hospital at Huazhong University of Science and Technology. All participants enrolled in this study signed written informed consent for participation, storage and use of surgically-removed cancerous and matched adjacent tissues. All methods of this study were carried out in accordance with relevant guidelines and regulations.

Data and tissue samples collection
We collected basic information such as age, gender, smoking status, smoking amount, and clinicopathological features including differentiated types of tumor Conclusion: IFT20 and GM130 protein expressions were negatively associated with tumor differentiated types, size, TNM stage and lymphatic metastasis of lung adenocarcinoma. Both IFT20 and GM130 proteins have some protective effects on the survival of lung adenocarcinoma patients with specific clinicopathological features. Keywords: Lung adenocarcinoma, IFT20, GM130, Clinicopathological features, Overall survival cells, tumor size, tumor node metastasis (TNM) stage, lymphatic or distant metastasis, cell proliferation index (Ki67), chemotherapy, and radiotherapy. Smoking amount (pack-years) for each smoker was calculated as packs of cigarettes per day multiplied by years of smoking. The differentiated types of lung adenocarcinoma were classified into micropapillary and solid, acinar and papillary, and lepidic types, which are identified as poorly, moderately and well differentiated tumors according to the criteria established by World Health Organization [17]. The TNM stage of patients with lung adenocarcinoma was diagnosed based on the 7th edition (2009) of the American Joint Committee on Cancer for Lung Cancer [18]. And the cancerous and matched adjacent lung tissues of each patient were collected for histological analysis.

Tissue chips preparation
Small pieces of tissue samples were fixed in 4% paraformaldehyde for 24 h, and then dehydrated, paraffin-embedded, sliced and stained with hematoxylin and eosin (HE), which were examined by microscope to determine the sampling location of the chip. After heating the paraffin block, the target samples taken out by a sampler were inserted into the prepared paraffin block receptor hole in sequence. The samples were merged at 60 ~ 65 °C for 20 min using tissue chip fusion instrument (BP0100, Biossci Company). The melted block was then put into a wax mold, embedded in paraffin and cut into 4 μm thick sections.

Immunohistochemistry (IHC)
The sections were deparaffinized, rehydrated, rinsed with distilled water, and repaired at high temperature and pressure for 1 ~ 2 min. After cooled to room temperature and washed with Tris-buffer solution (TBS) three times, the sections were incubated with fresh 3% hydrogen peroxide at room temperature for 20 min and blocked with 10% normal goat serum for 20 min. The sections were then incubated with primary antibodies (IFT20, 1:200, proteintech, 13,615-1-AP; GM130, 1:200, abcam, ab52649) overnight at 4 °C. After incubated at room temperature for 15 min, the sections were washed by TBS three times and then incubated with secondary antibodies (50 μl DAKO) for another 25 min at room temperature. Then the sections were washed three times, stained with diaminobenzidine and counterstained by hematoxylin.
The expressions of positive staining IFT20 and GM130 protein were indicated by the mean optical density (IOD/ area), the rate and intensity of positive staining cells. The average IOD/area and the rate of positively staining cells in five representative views (original objective 400×) from each section were analyzed by Image J software version 1.2.4 (the National Institutes of Health, NIH free software, Bethesda, MD, USA). The average IOD/area of these views represented the relative expressions of positive staining IFT20 and GM130 protein, and the rate of positive stained cells was equal to positive cells number/ total cells number× 100%. The staining intensity of positive cells of each section was denoted by staining intensity score (0 = no staining; 1 = weak staining; 2 = intermediate staining; 3 = strong staining), which was independently evaluated by five professionally trained persons without knowledge of the clinicopathological data of patients. The final score was determined as the consistent results made by three or more persons. D value is equal to the difference value of protein expression between cancerous tissue and matched adjacent tissue. The patients were divided into negative and positive groups based on the D values of IFT20 or GM130 protein staining intensity score, D value ≤0 was defined as the negative group, whereas D value > 0 was regarded as the positive group.

Statistical analysis
Continuous variables were divided into normal and nonnormal distribution, which were compared by Student's t test and Wilcoxon rank sum test, respectively. Categorical variables were analyzed by chi-square test or Fisher's exact test. The correlation between IFT20 and GM130 protein was assessed by Spearman's rank correlation. Associations of IFT20 and GM130 protein expressions with clinicopathological features of lung adenocarcinoma patients were assessed by multivariate logistic regression models, adjusting for multiple potential confounders including age, gender, smoking status, smoking amount. The survival analysis of patients with lung adenocarcinoma after lung cancer surgery was performed by Cox proportional hazard regression models. All statistical analyses were carried out with SPSS version 26.0 (SPSS Inc., Chicago, IL), and P < 0.05 was considered statistically significant.

The basic characteristics of all participants
As all the patients were divided into two groups based on the negative and positive expressions of IFT20 and GM130 protein, the number and proportion of IFT20 and GM130 positive expression patients were 168 (71.5%) and 151 (64.3%), respectively, In Table 1, we found that the amount of smoking was significantly less in IFT20 and GM130 positive expression groups than that in the negative groups (P < 0.05). Compared with IFT20 negative group, the proportion of male, current smoking patients was lower in IFT20 positive expression group, which was statistically significant (P < 0.05). The percent of clinicopathological features such as T3 ~ T4 stages, N2 ~ N3 stages and lymphatic metastasis was also much lower whereas the average OS of patients was significantly longer in IFT20 positive expression group than those in IFT20 negative expression group (P < 0.05). Besides, the proportion of poorly differentiated type cells was significantly lower in GM130 positive expression group when compared with the negative expression group (P < 0.05). And the average size of tumor (4.32 ± 1.89 cm) in GM130 positive expression group is also considerably smaller than that (5.02 ± 2.67 cm) in the negative expression group (P < 0.05).

Expressions of IFT20 and GM130 protein in lung tissues of patients with lung adenocarcinoma
As shown in Fig. 1 A and B, the expressions of IFT20 and GM130 protein were both positively expressed in cancerous and matched adjacent lung tissues of patients with lung adenocarcinoma. By semi-quantitative analysis, IFT20 expression (IOD/area, rate of positive staining cells and staining intensity score) and GM130 protein staining intensity score was significantly higher in cancerous tissues than those of matched adjacent tissues (P < 0.001) ( Table 2). Additionally, we found that the IOD/Area and staining intensity score of IFT20 protein were positively correlated with GM130 protein corresponding expression in cancerous tissues (r = 0.223 and r = 0.492, P < 0.001) and adjacent tissues (r = 0.385 and r = 0.424, P < 0.001) (Supplemental Table 1).

Associations of IFT20 and GM130 protein expressions with clinicopathological features of patients with lung adenocarcinoma
To illustrate the relationships of IFT20 and GM130 protein expressions with clinicopathological features, multivariate logistic regression models were analyzed    by adjusting for multiple potential confounders such as age, gender, smoking status, and smoking amounts. As shown in Table 3, each one-point increase in the D value of IFT20 protein staining intensity score was significantly associated with a 32% and 29% reduced risk for TNM stage in II ~ IV and lymphatic metastasis of patients, respectively (P < 0.05). The adjusted odd ratio (OR) and 95% confidence interval (CI) for rate of IFT20 positive cells in the highest quartile was 0.44 (0.20, 0.93) in comparison with that in the lowest quartile for patients with lymphatic metastasis. Although there was no association between IFT20 protein with types of tumor cell differentiation, tumor size or distant metastasis, each one-point increase in the D value of GM130 staining intensity score significantly reduced the risk of poorly differentiated type and tumors size> 7 cm by 29% and 38%, respectively (P < 0.05). In categorical analysis, there was also statistically negative relationship of GM130 protein expression (IOD/area or rate of positive cells) with the risk of poorly differentiated type cells, tumor size > 7 cm or TNM stage in II ~ IV. Furthermore, elevated rate of GM130 positive cells was also significantly associated with decreased risk of poorly differentiated type cells in a dose-dependent manner (P trend < 0.05). However, no significant relationships were observed between GM130 expression (IOD/ area, rate of positive cells or staining intensity score) and clinicopathological features of lymphatic or distant metastasis (P > 0.05).

Influence of IFT20 and GM130 protein expressions on the survival of patients with lung adenocarcinoma
Cox proportional hazard regression models were conducted to evaluate the influence of IFT20 and GM130 protein expressions on the survival of lung adenocarcinoma patients after surgery. As shown in Fig. 2 . 3 and Supplemental Fig. 1). In stratified analysis, enhanced IFT20 staining intensity score was significantly decreased the risk of death by 16% for patients without distant metastasis. And elevated the IOD/area of GM130 expression significantly decreased the death risk of lung adenocarcinoma patients with tumor size > 7 cm or distant metastasis by 54% and 65%, respectively (Table 4).

Discussion
In this study, we found the proportions of some clinicopathological features including T3 ~ T4 stage, N2 ~ N3 stage and lymphatic metastasis were significantly lower in patients with IFT20 positive expression than those with negative expression. And the average OS of patients was much longer in IFT20 positive expression group. Both the proportion of poorly differentiated type cells and the tumor size were significantly less in GM130 positive expression group when compared with the negative expression group. Multivariate logistic regression analysis was further confirmed that increased IFT20 protein expression was significantly associated with TNM stage in II ~ IV and lymphatic metastasis, whereas there was a significantly negative relationship between GM130 protein expression and poorly differentiated type cells and tumors size> 7 cm. And Cox proportional hazard models exhibited that poorly differentiated type cells, tumor size > 7 cm, TNM stage II ~ IV, lymphatic and distant metastasis were the risk factors for the survival of patients with lung adenocarcinoma. Consistent with our findings, a large amount of evidence has revealed that patients with lung adenocarcinoma have a poor prognosis when the tumor was poorly differentiated, in advanced TNM stages and with distant metastasis [19][20][21][22]. The micropapillary and solid subtypes classified as poor differentiation seems to be the most adverse prognostic factor [19,20]. The IFT20 and GM130 protein expressions were only negatively associated with the death risk of patients with specific clinicopathological feature such as tumor size > 7 cm or distant metastasis. These findings could indicate that (See figure on next page.) Fig. 2 The survival curves of clinicopathological features and IFT20 and GM130 protein expressions for patients with lung adenocarcinoma. *P < 0.05, compared with the blue line which is defined as reference line. The only one blue curve means that the survival curves for the Q2, Q3 and Q4 levels of IFT20 and GM130 protein overlapped with the reference line, respectively, after adjusting potential confouders such as gender, age, smoking status and smoking amount IFT20 and GM130 protein expressions may indirectly affect survival time by influencing some clinicopathological features which are associated with the survival of patients with lung adenocarcinoma. Primary cilium, a sensory appendage on the cell surface has a strong role in cell cycle control, differentiation and polarity that is correlated with prognosis of adenocarcinoma [23]. IFT20 is the smallest subunit of IFT complex, which is crucial for the formation of primary cilium within a cell [8]. Our results showed that the expression of IFT20 protein was significantly higher in cancerous tissues than that in matched adjacent tissues. A study reported that primary cilia frequency was significantly elevated in lung adenocarcinoma when compared with normal tissues [6], which may be partially caused by increased IFT20 expression in cancer cells. However, we found a significantly negative association of IFT20 staining intensity score with clinicopathological features like TNM stage and lymphatic metastasis of patients with lung adenocarcinoma. It is suggested that IFT20 might be a promoter of tumorigenesis, whereas it functions as a negative regulator in the progression of invasion and metastasis of lung adenocarcinoma cells. These seemingly paradoxical effects of IFT20 could be explained by a dual role of primary cilia mediated Hh signaling activation and suppression in oncogenesis and tumor progression [24]. Consistent with our results, Yang et al. reported that the level of IFT20 protein was negatively correlated with the malignancy of breast cancer cells, in which deletion of IFT20 promoted epithelial-mesenchymal transitions (EMT) and enhanced the migration of breast cancer cells by decreasing F-actin associated protein Transgelin-2 (TAGLN2) expression [25]. Overexpression of TAGLN2 could suppress lung metastasis in a mouse model of breast cancer [26]. A previous publication has identified that TAGLN2 was upregulated in patients with lung adenocarcinoma, most of whom were in the early stage of lung adenocarcinoma [27]. It is thus possible that negative relationship of IFT20 expression with TNM stage and lymphatic metastasis due to its role in regulating elevated synthesis of TAGLN2 in the early step of carcinogenesis, a hypothesis that will need to be explored in detail in the future.
A recent study has proposed that invasion of lung adenocarcinoma is associated with loss of cell polarization [28], which is known as a typical feature of EMT that promotes epithelial tumorigenesis and metastasis [29]. Some studies have indicated that the Golgi complex has been found to play a key role in promoting cell polarization [13,14]. As a Golgi marker, GM130 protein is involved in the maintenance of Golgi structure and the process of intracellular proteins transport to regulate the orientation of cell migration and polarity [15,16]. Depletion of GM130 resulted in a reduction of E-cadherin expression in epithelial cells, which  is considered to be an indicator of loss in cell polarity and epithelial identity [15]. It has been confirmed that a reduction in membranous E-cadherin/catenin complex expression was strongly associated with poor differentiation, distant metastasis and post-operative recurrence of patient with lung adenocarcinoma [30,31]. Therefore, the role of GM130 in regulating E-cadherin expression of epithelial cells may be used to explain our results of negative association between GM130 expression and the risk of poor differentiation, tumor size and shortened OS of patients with lung adenocarcinoma, which would be interesting to study further.
To the best of our knowledge, this is the first study to investigate the associations of IFT20 and GM130 protein expressions with clinicopathological features and survival of patients with lung adenocarcinoma. However, there are two limitations. Firstly, this is a case study in which the number of subjects is relatively small. But each patient's matched adjacent tissue could be considered as a self-control, which can reduce the bias to some extent. Secondly, although both IFT20 and GM130 proteins are localized in the Golgi apparatus and interact with each other, their interaction and mechanism on the development of lung adenocarcinoma could not be evaluated. We only found that IFT20 protein was associated with TNM stage in II ~ IV and lymphatic metastasis, whereas GM130 protein was correlated with poorly differentiated type cells and tumors size > 7 cm. The influences of these two proteins interaction on clinicopathological features of lung adenocarcinoma need to be explored in animal experimental models or case-control studies in the future.

Conclusions
In conclusion, we found the expressions of IFT20 and GM130 protein were negatively associated with clinicopathological features including tumor differentiated types, size, TNM stage and lymphatic metastasis of lung adenocarcinoma. IFT20 and GM130 proteins have some protective effects on the survival of lung adenocarcinoma patients with specific clinicopathological features such as large tumor size and distant metastasis. Thus, it is indicated that IFT20 and GM130 protein could play crucial roles in the development of lung adenocarcinoma.

Supplementary Information
The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12885-022-09905-6. Table 1. The correlation between IFT20 and GM130 protein in cancerous and adjacent tissues. Supplemental Fig. 1. The forest plot of Cox proportional hazard regression models with TNM stage included. To exclude collinearity, the three types' expressions of IFT20 and GM130 protein (IOD/area, rate of positive cells, and staining intensity score) were separately included in the model.

Learn more biomedcentral.com/submissions
Ready to submit your research Ready to submit your research ? Choose BMC and benefit from: ? Choose BMC and benefit from: